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Introduction 

The design of RF pulse sequences to produce selective excitation is a problem of prin- 
cipal importance in ajl applications of nuclear magnetic resonance. The oldest method 
of pulse synthesis is the Fourier transform method. Notwithstanding the fact that this 
method is an approximation for all non-zero flip angles, it gives surprising good results 
for flip angles up to about J. Over the past fi fteen years two, more exact approaches to 
solving this problem have been introduced and explored The slightly older, generally 
accepted method is the Skinnar-Le Roux or SLR algorithm introduced independently 
by M. ShinnaT and his co-workers and P. Le Roux. See for example [15], [16], or [12]. 
A more complete list of the original papers can be found in the references to [10]. 
The second method is to use the inverse scattering transform or 1ST which was first 
introduced as a technique for pulse synthesis in [8] and later refi ned in [1 1, 5, 6]. 

From the practical perspective either method of pulse synthesis is rather demanding 
and requires considerable computation. From the theoretical perspective the problem 
of pulse synthesis is very easy. Indeed, it is clear from the inverse scattering approach 
that there is an infinite dimensional space of RF-envelopes that will produce arty desired 
magnetization. The "real" problem of pulse synthesis is therefore to fi nd an "optimal" 
RF-envelope to produce a given magnetization profi le. In this paper we explore the 1ST 
approach and show how it can be used to design RF-envelopes which utilize the mini- 
mum possible energy for a given magnetization profi le. At the center of our discussion 
is a formula expressing the energy in the RF-envelope explicitly in terms of the mag- 
netization profi le and certain/ha? parameters that arise in the 1ST method. Formula? of 
this type were discovered, m the context of inverse scattering for the Korteweg de Vries 
equation, by Zakharov and Faddeev in 1971 ! See [19]. The formula we use appears in 
a 1974 paper of Zakharov and Manakov, see [20]. This formula is useful both for the 
design of pulse sequences as well as for their analysis, as it allows for an a posteriori 
determination of the erB ciency of a pulse sequence. 

Underlying the derivation of the SLR algorithm is the so called hard pulse approx- 
imation. One designs a pulse envelope that is a fi nite sum of ^-pulses which produces 
an approximation to the desired magnetization profile. As 5 -pulses are non-physical, 
requiring infi nite energy, an SLR pulse envelope is implemented by replacing the train 
of 5-pulses by a train of boxcar pulses. The 1ST algorithm instead constructs a smooth 
pulse envelope which again produces an approximation to the desired magnetization 
profi le. While such pulses are again implemented as a fi nite train of boxcar pulses, the 
approximations used in the 1ST are tax less singular than those used in the SLR. Indeed 
the 1ST provides more direct and greater control over the fi nal pulse envelope than 
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the SLR. When I undertook this work the central mystery forme was why Shinnar-Le 
Roux has remained the industry standard and the inverse scattering transform is largely 
ignored In addition to presenting the results, which seem interesting in and of them- 
selves, my second purpose in writing this paper is to present the JST approach in a 
more explicit and palatable form. This accounts, in part, for the somewhat expository 
nature of our presentation. 

In the fi rst and second sections we recall the connection between the Bloch equa- 
tion without relaxation and their spin domain formulation. In inverse scattering the spin 
domain Bloch equation is called the Zakharov-Shabbat 2 x 2-system or 2S-system. We 
next review the scattering theory for the ZS-system and relate the scattering data to the 
magnetization profi le. These sections closely follow [1 1], [10] and (I]. In the thiru" sec- 
tion we outline the inverse scattering transform via the Gel' fend-Levitan-Marchenko 
equation, or more briefly the Marchenko equation, following the treatment in [1], In 
the fourth section we state the formula relating the energy in the RF-envelope to the 
magnetization profi Ie and the location of the bound states. In the fifth section we dis- 
cuss the problem of implementing the 1ST. Under simplifying assumptions, we give 
a practical algorithm for solving die Marchenko equation. These assumptions are re- 
moved in Appendix A, where a modifi ed algorithm is given which applies in all cases. 
In the sixth section we address the problem of rephasing and the physical signifi cance 
of the time parameter. The seventh section contains several examples. Inthe eighth sec- 
tion we analyze the approximations used in both the SLR and 1ST approaches. Proofs 
of several mathematical results are contained in Appendix B. 



1 The spin domain Bloch equation and the problem of 
RF-puIse synthesis 

The Bloch equation without relaxation is usually written in the form 

dM 

— = yMxB. (1) 

Here M is the magnetization, B is the applied magnetic fi eld, / is time and y is the 
gyromagnetic ratio. Since a vector evolves with constant length under this equation, 
the fundamental solution £%(*) can be regarded as *S0(3)-vaIued function of time. 
Here we normalize so that Ub(0) is Id 3 , the 3 x 3 identity matrix. If M(Q) = Af 0 
then Af (0 = t/ 0 (f) Afo* The solution to the Bloch equation is linear in the initial data. 
Throughout this paper we assume mat solutions of the Bloch equation are normalized 
to have length equal to one. 
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The Bloch equation,* usually analyzed in a "rotating reference" fiame. Ordinarily 
the rotating reference frame is related to the "laboratory frame" by a time dependent 
orthogonal transformation of the form 



Fit) 

so that 



[cos0(f) -sin0(f) <fj 
» sin0(/) cos0(/) 0 J t 
L 0 0 lj 



mo = <F(o«<0. , (2) 

We use m to denote the magnetization in the rotating reference frame. Larmor's theo- 
rem implies that if M satisfj es (1) men m satisfi es 

dm 

— = yi* x Heff (3) 

where 

Beo(t) = F~ x (/)*(/) + with fi(f) = [0, 0, 0'(fl]t. (4) 

Rotational remark: Most of the vectors used in mis paper are to be thought of a col- 
umn vectors. The notation [a, b t c$ refers to the transpose of the row vector which is 
therefore a column vector, i.e. 



In most applications of this method the function 9(0 is selected to render the z- 
component of B e tr independent of time, 

BcffU; t) == (and), on(t) 9 y l f). (5) 

The constant value / is called the offset frequency or resonance offset If we set cb(t) = 
<o\ (0 + ianit) then, in the laboratory frame, the RF-envelope takes the form: 

Pit) = Ute&OWWK bn&(f)* m 9 0]t. 
The energy in the RF-envelope is given by 

op 00 / 

\P(f)\ 2 dt= f m)\ 2 dt. (6) 
-00 -00 

The magnetization profi le is a unit vector valued function defi ned for / € U f 
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In essentially all MR applications m°°(f) = [0,0, 1]* for / outside of a bounded 
interval. The problem of RF-pulse synthesis is to fi nd a time dependent complex pulse 
envelope c5(0 so that, if J*efK/) is given by (5), then the solution of 

dm 

t) = Y'»(f ;t)xB e{ t(f ;0 (7) 

with 

Jta >»(/;/) = [0.0, l]t (8) 

satisfi es 

fa£e-V'(mx + /m 2 )(/; t), m 3 (/; I)] = + /72~(/)]. (9) 

We have used the standard complex notation for the transverse components of the mag- 
netization. If a> is supported in the interval [/o. ' 1 ] theft these asymptotic conditions are 
replaced by 

«(/;*<>) = [0,0,1]*, 
le~V'* On, + im 2 )(f ; *i), w 3 (/; = [(mf + imf )(/), m§°(/)]. 

The mapping from <y to is highly non-linear. 

To solve the problem of KF-pulse synthesis it is convenient to introduce the spin 
domain formulation of the Bloch equation. Instead of a unit vector m in (R 3 we solve 
for a unit vector * in C 2 . This vector satisfi es the 2 x 2 matrix equation: 

i—ZQ-wt. (ID 

Here Q = -[y^i (/). y<B2(0» /] and <r are the Pauli spin matrices: 

— B i]- — p »]• —p. ..]• 

Assembling the pieces we see that * satisfi esf 
with 

* = 2 * 98S "2" W,) " - ^(O). 
A simple recipe takes a solution of (12) and produces a solution of (7). If *(f ; /) = 
(I; f). lte(£; 01* satisfi es (12) then the 3-vector 

«(/; 0 - BR«(#i*«, 2im(^fifr 2 ). 1*1 1 2 - m\ 2 tf(j; 0 (13) 



1 We follow the standard practice in the MR literature of using z* to denote the complex conjugate of the 
• rz. 
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satisfi es (7). If in addition 

then m satisfies (8). Thus the RF-pulse synthesis problem is easily translated into a 
inverse scattering problem for equation (12). This is described in the next section. 
Following the standard practice in inverse scattering we refer to equation (12) as the 
ZS-system and q as the potential. 

2 Scattering theory for the ZS-system 

Scattering theory for an equation like (12) relates the behavior of lim,-*^ ; /) to 
that of lim^oj ^(£ ; /). If q has bounded support then the functions 

m- [a] 

are a basis of solutions for (12) outside the support of q. If the L^nonn of q is fi nite 
then it is shown in [1] that (12) has solutions which are asymptotic to these solutions 
as/ -> rfcoo. 

Theorem 1. If \q ]| L i is finite then for every real f there are unique solutions 

\fr 2 +(£) and ^,_(£ ), ^ 2 _(£) 

to equation (12) which satisfy 

J^*'*l-«;<> = [i]. ^e-'^^O^f. 0 ,] (14) 

,!™*'<V I+ tt;0 = [j], < ita e -'^ 2+ «;/) = [;] (15) 

Thesolutions lh_(g), # 2 +(£) extend as analytic functions off to the upper half plane, 
Irn £ > 0 flnrf (f ), \fr 1+ (f ) extend as analytic/unctions off to the lower half plane, 
Im£ < 0. 

The proof of this theorem can be found in [1J. We briefly outline the construction 
of these solutions as we will need these formulae later on. 

Outline of proof The idea is to re-express the boundary value problems for the differ- 
ential equation as Volterra integral equations. Let 

'■-[&]• *--[&]• 



(16) 
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These solutions satisfy: 

— oo 
t 

-00 

/ 

lk*2-«:0*~'*' = -l + y A/2(?U,5)fe-($;5)e-^, 
-oo 

I 

*i2-e;0e" ,|# = - y e- 2/f( '-* ) g(5)^22-(?;5)e- , ^^ l 

-00 

where 

A/i(|; /,*) = -**(*) y ^-^(x)rfx, 
i 

M 2 tf; f,*) « -*(*) y f**-*>4*QQdx. 

a 

The solutions nonnalized at +oo satisfy similar equations. If we set 

t 

a©) - y if«k& 

— oo 

then it is shown in [1] that these equations can be solved for § in the appropriate half 
plane by iteration and 

Itli-G; ry*'l < /o(20,(O). Itea-fc; 0*-*| < AC2fii(0). (18) 

Here /o is the classical /o-Bessel function. The same argument applies mutatis mutan- 
dis to obtain $ 1+ and #2+* □ 

For real values off the solutions nonnalized at -co can be expressed in terms of 
the solutions nonnalized at -Hx> by the linear relations 

0 = *<?)*!+(*; 0 + *(*)*2+<S; 0. nm 

The functions a, & are called the scattering coefficients for the potential q. The 2x2- 
matrices [^,_^ 2 -] . [^1+^2+] satisfy 



<ft O °~f 3 T ,iS i 3 O 5? 
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The scattering matrix for the potential 9 is defi ned to be 

Recall that the Wronskian between two C 2 -valued functions is defi ned by 
W{U. »)(*) L H|(0V2(0 - u 2 (/)vi(0. 

If «(0 and i7(/) ate solutions of (12), for the same value off, then JF(«(0. „(,)) is 
independent of/. It is not diffl cult to show that ' 

«tt> = [*n-(f: Ofe»<-tt; 0 - tei-(f ; Oi^ia+tf ;/)] » IT(^,_. * 2+ )(0, 
*«) = -Wu-<£; Otoi+(f ; 0 - toi-(f ; 0*n+(f ; 01 = f l+ Xt). 

(22) 

It follows from (22) that a extends to the upper half plane as an analytic function. If a 
has an integrable derivative then, using (16) and (22) it can be shown that 

00 

0(M } = 1 + 27| / ,9(S)|2 * + °(^> (23) 
-oo 

as |f | tends to infinity in Imf > 0. As W{^_, f 2 _) = -lit follows that 

l«(f)l 2 + |A(f)l 2 = l. (24) 

and therefore (23) implies that 

WS>I = 0(~) as f ±oo. (25) 

These results are proved in (1 J. 

Assuming that t Jq{t) is integrable for aU j it is shown in [1] that a has finitely 

many zeros fata* > 0. tet (f, M be a list of the zero's of a. FoHach / ? 

follows from (22) that there is a non-zero complex number Cj so that 

*!-(£/) - Cjir 2+ %), y=i Nt (26) 

The constants {C,, . C N ) appearing in (26) are called norming constants. Equa- 
tion (12) can be rewritten fa the form ^ 

From tins formulation it is clear that f should be regarded as a spectral parameter If 
f has positive imaginary part then *,_(£ ; 0 decays exponentially as / tends to -co 
Ml' 'V^" 3 ™ exponentially as / tends to +oo. If ««,) = 0 then (26) implies 
««>at *,_(£,; 0 decays exponentially at both ±oo and therefore the function *,_(£• ») 
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2*7 y 



belongs to L 2 (K; C 2 ). Thus the operator on the left hand side of (27) has bound states 
for these values of the offset frequency. 

We generally assume that the zeros of a are simple and that their imaginary parts are 
positive. This is mostly,to simplify the exposition, there is no diffi culty, in principle, if 
a has real zeros or higher order zeros. Using the argument principle and (23) it follows 
that 

am • (28) 

Definition 1. The pair of functions (a(£), 6(f)), for £ e R and the collection of pairs 
KSjfCj) : ; = I,..., AT) define the scattering data for equation (12). 

The scattering data are not quite independent Since a is analytic in the upper half 
plane it follows from (24) and (23) that, for Im f > 0, we have: 

■«>-n(^)-[i/!«^]. „ 

see [7]. The reflection coeffi cient is defi ned by: 

A priori the reflection coeffi cient is only defi ned on the real axis. We rewrite (29) in 
terms of r : 

° a, =}3(^H^/#H- <3,> 

Both (29) and (31) have well defi ned limits as £ approaches the real axis. 

Definition 2. The function r(f), forf e Rand the collection of pairs {(£, , C, ) : j = 
1 ..... N) defi ne the reduced scattering data for equation (12). 

Implicitly the reduced scattering data is a function of the potential q. In inverse 
scattering theory the data {r; ($1,0,..., (£„, Cat)} are specfi ed and we seek a po- 
tential q which has this reduced scattering data. The map from this data to q is often 
called the Inverse Scattering Transform or 1ST. 

We close this section by rephrasing the RF-pulse synthesis problem as an inverse 
scattering problem. Recall mat the data for the pulse synthesis problem is the magne- 
tization profi le nf° which we now think of as a function off = f/2. Using (13), the 
solution to the ZS-system defines a solution to (7), satisfying (8). It follows 
from (15) and (19) that 

0 ~ [^Jji^] as / +co. (32) 
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Therefore 



where, as before, we use die complex notation for the transverse components of mi-. 
If also satisfi es (9) then it follows from (33) and (24) that 

*(f ) _ (wu-+ftBai-)g;»X- at «« 
a(£) «-oo l+ m3 ,_g;/) 

= K°+i»if)($) < 34 > 

l+mf(?) ' 

In light of this formula we sometimes refer to r as the stenographic magnetization 
profile. Ifq has support in the ray (-co, /j] then 

l+m 3 i-(£;/) ^ 

is independent of / for t > t\. 

As m°° is a unit vector valued function we see that the reflection coeffi cient r (§) 
uniquely determines and vice-versa. Thus die RF-pulse synthesis problem can 
be rephrased as fee following inverse scattering problem: 



Find a potential q for the ZS-system so that the reflec- 
tion coeffi cient r satisfi es (34) for all real f . 



Note that the pulse synthesis problem makes no reference to the data connected to the 
bound states, i.e. {(£,, C/)}. Indeed these are free parameters in the pulse synthesis 
problem, making the problem higly underdetermined. Perhaps the main result of this 
paper is that, for a given magnetization profi le, the RF-envelope requiring the minimum 
energy is the one for which the ZS-system has no bound states. 
Remark 1. The material in this section is largely adapted from [1], [1 1] and [10]. 



3 Inverse scattering for the ZS-equation 

There are several different approaches to solving the inverse scattering problem stated 
at the end of the previous section, see [1], [2], [4] or [7J. In this section we describe the 
solution of this problem via the Marchenko equations. We again follow the presentation 
in [1]. We generally suppose that the {£, } are distinct complex number with positive 
imaginary parts. To begin with we assume that the reflection coefficient is smooth 
and rapidly vanishing at ±oo. Using limiting arguments the tatter restriction is easily 
relaxed. 

The Marchenko equations are systems of integral equations of Fredholm type 
Given the scattering data {r; (f ! , C,) C N )} defi ne the function 

/CO = ^T / rQ*"4 - i Cje*>'. (36) 
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This is the inverse Fourier transform of r with a correction added to account for the 
bound states. Note that the fi nite sum is exponentially decreasing as / tends to +00. 
Defi ne the 2 x 2 matrix valued function: 

*»=U> 

For each number t <= R we seek a 2 x 2-matrix valued function defi ned for 

5 € [r, 00), which solves the integral equation: 
00 

3£<&) + y %(<r)S?(y + <r)</<r » -g?(f + j?) for* e [r, 00). (37) 

Under our assumptions on r it is not diffi cult to show that, for every t e U, (37) is a 
Fredholm integral equation on Z 2 ([r, 00) ; o/ 2 ). Since 9{t+ .) belongs to Z 2 ([r, 00); gfe), 
the solvability of (37) follows if it can be shown that any solution to the homogeneous 
adjoint equation, 

00 

£(s) + j <£(<x)&>(s + a)dv = 0 

i 

is identically zero. This proof of this statement can be found in [1J. The uniqueness 
of the solution to (37) follows from the fact that homogeneous equation corresponding 
to (37) also has only the trivial solution. 

Remark 2. For the theory of Fredholm integral equations see [18], For the moment we 
note that the known proofs for the existence of a solution to a Fredholm equation do 
not generally provide a practicable algorithm to fi nd the solution. 

The connection between the inverse scattering problem and the Marchenko equa- 
tion is summarized in the following theorem: 

Theorem 2. Given a smooth, rapidly decaying reflection coefficient r and a finite set 
of pairs {(£,, Cj) : j e> l f . . . , JV), with the distinct, lmfy > 0 and Cj £ 0 for 
j = 1 , . . . , Jv\ the equation (37) has a unique solution for every t GU.If 

*(0=-29C,*(0 (38) 

then the ZS-system with this potential has reflection coefficient r. It has exactly N bound 
states for frequencies and the relations (26) hold at these points. 

This theorem is proved in [7] and, in a slightly different formulation in [2]. 
Remark 3. As noted above, there is a similar result if a has higher order zeros or zeros 
on the real line. In this case the defi nition off needs to be modifi ed, we return to this 
in Appendix A. 

We defer a discussion of the practicalities of solving (37) to section 5. For the 
moment suffi ce it to say that this equation can be solved by a simple iteration for values 
off such that 

00 

j mender < I, 
21 
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while a somewhat indirect method is required for smaller values of /. There is also 
"left" Marchenko equation, given in section 5, where the integrals are over half lines 
of the form (-oo, r]. In this formulation the fi nite sum of exponentials is decreasing as 
/ tends to -co. 

The regularity and decay of 9 are related to the regularity and decay of r exactly 
as if they were a Fourier transform pair. The following result is a special case of a 
theorem proved in [2]: 

Theorem 3 (Beals-Coifman). Suppose that r is a reflection coefficient and q is the 
corresponding potential for the ZS-system. Let k be a positive integer. 



then the distributional derivatives dfqfor 0 < j < k belong to £ 2 (R) + Z,°°([R). 
2. If the distribution derivatives d{r e L 2 (R)for 0 < j < k + I then t k +*q e 

2, 2 (R) + L°°my 

Remark 4. There are no smallness hypotheses in this theorem. Note also that this theo- 
rem does not mention the bound states. While the bound states exercise a strong effect 
on the detailed structure of q, their presaence or absence has no effect on the gross 
smoothness and decay properties of the potential. These features of q are governed 
entirely by the reflection coeffi cient. 

Remark 5. It has been noted many times in both the inverse scattering and MR litera- 
tures that, as q tends to zero, r*(f ) tends to £(-2|). The results of Beals and Coifman 
shows that the connection between the inverse scattering transform and the Fourier 
transform is much deeper than that previously described. There is also an analogue of 
the Paley-Wiener theorem (see Appendix B): If the reflection coefficient has bounded 
support then the potential does not This is easier to see in the logically equivalent 
contrapositive form: If the potential has bounded support then the reflection coeffi cient 
does not When the potential has bounded support the scattering coeffi cients (a, b) ex- 
tend, as analytic functions, to the entire complex plane. This means that the set where 
b and therefore r vanishes is discrete. Hence the reflection coeffi cient, restricted to the 
real axis, cannot have bounded support 

Remark 6. Beals and Coifman show that the 1ST can be extended to square integrable 
reflection coeffi cients such that %r also belongs to L 2 (K). By using smooth, compactly 
support approximating sequences, the estimates in the theorem can be extended to data 
with fi nitely many derivatives and fi nite rates of decay. 

The results of this section show that, at least in principle, the 1ST provides many 
solutions to the RF-pulse synthesis problem. A practical magnetization profile has 
bounded support and may be approximated by a smooth function. This leaves two 
essential problems: 
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1. Given (and therefore r) how should die bound states be selected to give 
an "optimal" RF-pulse envelope which will produce the desired magnetization 
profile? 

2. How is the Marchenko equation to be solved? 

The notion of optimality depends on the intended application. It usually entails a 
balance among the energy, duration, oscillation and need for rephasing of the RF- 
envelope. In the remainder of the paper we focus on the construction of the minimum 
energy RF-envelope. 

4 The energy of the RF-envelope 

We now state a formula for the energy of the pulse envelope in terms of the scattering 
data. The underlying results are due to Zakharov, Faddeev and Manakov. 

Theorem 4. Suppose that q is a sufficiently rapidly decaying potential for the ZS- 
system with reflection coefficient r and discrete data{% % Cj\ j = 1 N} then 

oo oo 

\q{t)\ 2 dt = - J log(l + |rG)|V* + 4£lhlfj. (39) 
—op — oo y=i 

The proof of this result can be found in [1] or [7]. Note that the norrning constants 
play no role in this formula. Formula (39) is just one from an infi nite sequence of 
formulas relating functional of the potential to functionate of the reduced scattering 

Combining (6), (39) with (34) we obtain the following simple corollary. 

Corollary 1. If m°° is a sufficiently smooth magnetization profile such thatm?°+im?> 
vanishes outside a finite interval then the total energy of any RF-envelope co which 
produces this magnetization profile satisfies the estimate 

Equality holds in this estimate if and only if the ZS-system with the corresponding 
potential has no bound states. 

From the corollary it is evident that the lowest energy RF-envelope is obtained by 
solving the Marchenko equation with 

oo 
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Remark 7. There are several other immediate conclusions which can be drawn from 
this formula. Often times one seeks to specify a particular dependence for the phase 
of the transverse magnetization. Changing the phase of the transverse magnetization 
amounts to replacing (/»~(f ) + imf « )) by </*f <£) + im?(&)e**G\ where * is a 
real valued function. From (40) it is apparent mat the energy required for the minimum 
energy RF-envelope is independent of the phase of the transverse magnetization. Note 
however that the minimum energy pulse will in general require rephasing. This is 
addressed in section 6. 

Remark 8. Very important types of pulses in MR imaging are inversion or refocussing 
pulses. Such a pulse carries [0, 0, 1] to [0, 0. -1]. If this is the case at offset frequency 
£ then r(£) = co. It would therefore require infinite energy to exactly invert [0, 0, 1] 
for offset frequencies belonging to an interval of positive length. The energy required 
to carry [0, 0 t 1] to a vector of the form [m u m 2t -1 + e] for 0 < e < € and offset 
frequencies belonging to a band of width B f is at least 



Ideally we would like to have a functional that estimates the support of q in terms 
of the scattering data. An explicit formula of this sort has proved elusive. On the other 
hand, as we shall see in section 5, the Marchenko equation implies that the effective 
support of q is largely determined by the effective support of the function / defi ned 



In the absence of bound states it is not hard to show that if r(£) is the reflection 
coeffi cient determined by q(t) then r(A£ ) comes from the potential Xr l qOr*Q. This 
agrees with the well known heuristic principle that the effective support of the RF- 
envelope is inversely proportional to the bandwidth of excitation. It also shows that the 
maximum amplitude and energy of the RF-envelope are proportional to the bandwidth. 
The L -oscillation of the RF-envelope, / \q t \ 2 dt t is proportional to the square of the 
bandwidth. 



5 Solving the Marchenko equation with no bound states 

In this section we discuss an algorithmic approach for obtaining the minimum energy 
RF-enyelope. It entails solving the Marchenko equation, (37) and a related equation, 
(42) given below. Our method involves solving the Marchenko equations directly via 
iteration. As noted aboye, this requires some care as large data may lead to a divergent 
series. In [11] Rourke and Morris introduced a method for solving the Marchenko 
equation by adapting a method of Moses and Proesser, see [9J. The analysis in [11] 
differs in one small particular from that in [9], leading to solutions of the Marchenko 
equation with non-trivial bound states and therefore sub-optimal energy requirement 

The method of Moses and Proesser is to approximate the reflection coeffi cient by 
a rational function with all of it poles in the lower half plane. This ensures that the 
resultant potential is supported in (-oo, 0] and has no bound states. In [1 1] this condi- 
tion on the placement of the poles is not respected Indeed in examples 2 and 3 of [1 1] 




in (36). 
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they use a rational approximation to r with poles in both half planes; the solutions they 
found were supported in a proper subinterval of (-00. 0). However the poles of r in 
Ae upper half plane force the resultant potential for the ZS-systera to have non-trivial 
bound states. These increase the needed energy and oscillation without affecting the 
magnetization profi le. Regrettably the method of Moses and Proesser is not generally 
applicable to fi nd the minimum energy RF-enveippe. This is due to the well known con- 
sequence of the maximum principle in the theory of functions of a complex variable: 
A function on the real line which tends to zero at infi nity is only well approximated 
by rational functions with poles in the lower half plane if it is the boundary value of 
an analytic function defi ned in the upper half plane. In light of a classical theorem of 
Fatou, this excludes any function supported in a bounded interval, see [13], 

We are left with little choice but to solve the Marchenko equation directly. We 
now present an algorithm for doing this. The algorithm is described in the realm of 
continuum mathematics but can easily be approximately implemented on a computer. 
We will return to the question of implementation in a later publication. After stating the 
"left** Marchenko equation we give an estimate which shows that, for a given reflection 
coeffi cient, the Marchenko equations can be solved for suffi cienty large |/|. There are 



1 . Using the Marchenko equation we solve for the potential one tat a time. 

2. If q = qi + q2 where q\ and q 2 have disjoint supports then the scattering data 
for q has a simple, explicit expression in terms of the scattering data for q\ and 

Using these ideas leads to a method which allows us to work inward from ±00 to 
determine the potential for all values oft. 

Theorem 2 and its extension to less regular data imply mat, given a sufficiently 
regular reflection coeffi cient r, there exists a potential q so that the ZS-system, with 
this potential, has reflection coeffi cient r and no bound states. In this and the following 
section the data r and the corresponding potential q are fi xed. We assume that r decays 
rapidly enough and is suffi ciently regular for q to be bounded and integrable with two 
continuous derivatives, see Theorem 3. 

Equation (37) involves integration over positive half rays. Because of the diver- 
gences indicated above it is useful to work both ends against the middle. For that 
purpose we give the left Marchenko equation which involves integration over negative 
half rays. For later applications we give the left Marchenko equation with fi nitely many 
bound states. Recall that r = b/a; given r we can use (3 1) to determine a and therefore 
0. Deri ne 

and let {f j , . . . , $ N ] be points in the upper half plane. There are non-zero constants 
[Cj } so that the kernel function for the left Marchenko equation is defi ned by 

8(t) = h f - 'E £ /' 9 ' (41) 

-00 7=1 



r 
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The exponential correction terms are exponentially decaying as t tends to — oo. The 
kernel for the left Marchenko equation is 



-**(0 oj' 
For each fixed* € K the equation is 



y ^/W < 5(x+5)^ = -^a + ^)forjr€(-oo f f]. 
— oo 



(42) 



If S€ t (s) solves this equation then the potential defi ned by 

q(t) = 2X t (t) 

has reflection coefficient r and bound states at {fy). In the remainder of our analysis 
we concentrate on the "right" Marchenko equation, (37), with the understanding that 
everything said applies mutatis mutandis to the left formulation as well. 

Let 1-U denote an operator norm on 2 x 2-matrices; this means that ffl-J is a norm 
which satisfi es: 

usi < mi m 

For example we can take 
Defi ne the functions 



Ml.nnJ^ 
v&> Dv|| 



M/(0 = supl9?(j)l, 
*>/ 
00 

'/»- f WisMds. 



(43) 



We say that the Marchenko equation can be solved by iteration on the interval It, oo) 
if the sequence defi ned by 

K?fr) = -SFCs + 0. 

. OO 

3f/ (5) . -&<s + 0 + j %{-\s + x)®(t+x)dx % j - 1, 2. . . . (44) 
r 

converges uniformly to 9£/(j) for jr € [f, 00). 

Theorem 5. #7/(2/) < 1 /ten the Marchenko equation can be solved by iteration on 
the interval [r, 00). The solution satisfies the estimate 
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The proof is given in Appendix B. Note that no assumption is made here about the 
presence or absence of bound states. 

The estimate (45) implies that, so long as the hypotheses of die theorem hold, the 
potential q satisfi es the estimate 



WV "~ l-//(20 



(46) 



This in turn explains our claim that the essential support of q is determined by the 
essential support of /. A similar result holds for the left Marchenko equation. Taken 
together they show that if r decays rapidly enough and is suffi ciently smooth then the 
potential q can be determined, by iteration outside of a bounded interval [/_, f + J. If 
//(-co) < 2 then q can be determined on the whole real line by iteration. 

If //(-co) > 2 then an additional trick may be needed. If we rewrite the potential 
q as a sum of potentials with disjoint support: 



q = X(-oo./ 0 W + X(/ 0 .ooj? 

= q\ + qz 



(47) 



then the scattering matrix of q is easily expressed in terms of the scattering matrices of 
qx and q%. Indeed we have two formulae. 

Theorem d. Suppose that q = q\ + q 2 where the support ofq\ is to the left of the 
support of q 2 . If 



are the scattering matrices defined byq i% i = 1,2, respectively then the scattering 
matrix defined by q is given by / 



a 2 -b* 
a 2 b* 



(48) 



These formula are elementary consequences of the defi nition of the scattering matrix. 

The abstract existence theory for the Marchenko equation assures that, given a 
sufficiently smooth, square integrable reflection coefficient r, there are solutions of 
the corresponding Marchenko equations for every t € U and therefore a potential q 
which defi nes a ZS-system with no bound states and the given reflection coeflfi cient In 
light of Theorem 5 there are numbers {/i 0) , tf>) so that we can solve the left and right 
Marchenko equations by iteration to fi nd 



(0) 

=Xl,f.oo)9- 



r «o,g and 



(49) 
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If »i 0) > »j 0) then wc are done. Otherwise we need to solve the forward scat- 
tering problem for the potentials a£°> and q™ to obtain the scattering coefficients 
{(a ± , 6J.0). Using Theorem 6 twice we obtain the relation: 

L* -«-.J-[*f arjK JL-^ «H (50) 

Here ( aj . 6,) are the scattering coeffi cients for the (unknown part of the) potential 

9i - X^m (51) 

Using (50) we solve for the scattering coeffi cients of qi . Since q x is a potential with 
bounded support the functions (a, , b x ) extend as analytic functions to the entire com- 
plex plane. It is not clear whether or not the ZS-system with potential a, has any 
bound states. In fact our algorithm applies in either case. To simplify the exposition 
we assume, for now, that for all a <fi the potentials Xla . P] q have no bound states, 
fa Appendix A we explain how to remove this hypothesis. The iterative step in the 
algorithm relies on the following Paley-Wiener type theorem: 

Theorem?. Supposethatq is a potential/or the ZS-system with support in the half line 

J*H!t ** e ^'fi"^" ffi>r the Marchenko equation, defined 
w (30), (g defined in (41)) issupportedin die ray (-oo, 2t+] ([2r_, oo)). 

The proof of this result is given in Appendix A. 

Assuming that q x has no bound states, the pair of functions (ai , 6,) is the scattering 
date for a potential with support in the fi nite interval (£ 0 \ ,<?>]. Theorem 7 therefore 
applies to snow that 

00 

/>(') = i/nCf^'rf* (52) 

— OO 

fa ajfonction with support in the interval (-oo, 2/<°>], with an analogous conclusion 
From Theorem 5 it follows that there exists another pair (fi 0 , f£ } ) such that 
/i 0> <ri'\ /i»<rf 

and the Marchenko equations, with the kernels defi ned by a , fx, can be solved by 

iteration on the intevals (-oo. ri«>l and [f<!\ oo) respectively. Note that the potential 

q is now determined on the complement of «£>]. Next the scattering data for the 
potentials 

9- } = X[,<o> , «)ja and 

„(0) " " (53) 
*+ = x t«i'Vf 1^- 

are determined. Using Theorem 6 this data can be used to determine (a 2 , A,) the 
scattering data for q 2 - X fl O) ,(i),a. ' 
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The method is applied tterativery to obtain sequences 

ri 0) < < £+0 < ... < ,(,+» < ,0) < ^ ^ 

In the Appendix A we show that, given suffi ciently smooth initial data r there is A > 0 
so that, for all j > 0, 

|^>«4f+»|>A. 

Hence after fi nitely many steps the potential q is determined on the whole real line. 

We close this section by stating our algorithm as a sequence of steps. In the next 
section we indicate the modifi cations needed to handle the possible occurence of bound 
states for the "intermediate" potentials [qj) t 

Step 0: 

Solve the left and right Marchenko equations by iteration on "maximal' 9 in- 
tervals (-co, /i°>], [/<?>. oo). If ri°> > we are done. Otherwise deterinine 
the scattering data for ?i 0) and q£ y . Apply Theorem 6 to determine the scat- 
tering data for q\ . Set j = 1 and proceed to step 1 . 

STEPl: 

Use the reflection coeffi cient for ^ in (52) to determine the kernel function 
fj (gj) supported on (-oo, 2t^" l) ] ([2/^~ ,) , oo)). 

Step 2: 

Solve the Marchenko equation on the interval [r£\ t%~ n ] (R?" 1 *, t^J). If 
> we are done, otherwise go to step 3. 

STEP 3: 

Determine the scattering data, on the real line for the potentials 
and apply Theorem 6 to determine the scattering data, on the real line, for 
Step 4: Go to step 1 with j replaced by j+1 . 

6 What is time? x 

Before we turn to the computation of examples we need to consider the physical sig- 
nificance of the f -parameter which appears in (1 1). This is intimately connected to the 
problem of rephasing. The main clue comes from (35). Suppose that q is a potential 
supported in ft), t\] with reflection coeffi cient r. After rime = 4 a solution to (1 1) is 
freely processing in the local (gradient offiset) ^o-field. 
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If /i — 0 then, at the end of the RF-excitation, the magnetization satisfi es 

(«ii-+/W2i-)tt;0) 

l+«3i-tt;0) 55 r(?) - < 54 ) 
We have therefore achieved the desired magnetization profile without any need for 
rephasing. If on the other hand *k 0 then, at the end of the excitation, the magnetiza- 
tion satisfi es 

In order to get the desired magnetization profi le the spins must freely process for (4 1- 
units of time. Finally if t , > 0 then the spins must process for -units of time. In 
other words either a 1 80°-refocussing pulse must be applied or the gradient needs to be 
reversed, followed by t\ -units of free precesion, in order to achieve the magnetization 
profi le specifed by r. 

In RF-pulse synthesis the data is specifi ed in the frequency domain. Given r (£) and 
perhaps some bound states {(£, , C,)) the 1ST produces a potential which is a function 
of /. The support of this potential is determined by the scattering data and thereby 
determines what sort of rephasing is needed to achieve the specifi ed magnetization 
profi le. If the essential support of the potential is contained in [©, *i] then, according to 
the sign of fi, me necessary rephasing is determined as in the previous paragraph. The 
main point to note is that the time dependence of the potential is an output of the 1ST. 
In particular, the numerical value of t u in appropriate units, determines the rephasing 
which is needed to achieve the specifi ed magnetization profi le. When computing the 
total energy requirement for an RF-envelope with t\ > 0 one must include the energy 
of the (hard) 1 80°-pulse or that of transients produced by reversing the gradient. 

Two tools are available to manipulate the support of q. On the one hand we can 
use smooth approximations to r so that its Fourier transform is rapidly decreasing. The 
estimate (46) shows that this reduces the essential support of q. The other tool is the 
introduction of bound states, coupled with methods of meromorphic approximation. 
Rourke and Morris approximate r by rational functions which decay at infinity and 
obtain potentials with support in (-oo. 0]. The rational approximations to r which 
they use have poles in the upper half plane and these in turn lead to an increase in the 
energy and the oscillation of the RF-envelope. This is the price for an envelope that 
does not require rephasing. For the standard pulses used in MR imaging the minimum 
energy pulse in general requires rephasing. 

The problem of finding approximations that carefully balance the rephasing re- 
quirements with the energy of the pulse itself seems a very interesting fi eld for further 
investigation. Through the Marchenko equations and the Paley-Wiener Theorem II 
it translates to the problem of approximating the function r, defi ned on real axis, by 
meromorphic functions p, defi ned in the upper half plane, which satisfy an estimate of 
the form 

oo 

/ e-*i T \p$ + it))\*d$ < C for j, > 0. 

— OO 

Here C and T are positive constants. The resultant potential will be supported in 
("~°°» n. We will return to this variational problem in a subsequent publication. 
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7 Examples 

8 Analysis of the hard pulse approximation 

The Shinnar-Le Roux algorithm uses the hard pulse approximation to design an RF- 
envelope. Given a magnetization profi le r <£) the SLR-algorithra designs a potential of 
the form 

N 

go = ~j A), 

to produce a good approximation to the desired profi le over a sub-interval of 

The reflection coeffi cient, jo defi ned by g> is aperiodic function of period A~ r 2jr, The 
SLR-algorithm provides an iterative method for detennining the coeffi cients [p,j }. 

Of course a sum of 5-pulses is non-physical, requiring infinite energy to realize 
The RF-envelope which is actually used is a "softened" version of q 0 . For example one 
could replace each pj&(t - jA) by a boxcar pulse of width A with the same area, 
leading to the softened pulse: 

TV 



7=1 



see [14]. There is no reasonable sense in which the difference q 0 - q\ is "smalL" 
Nonetheless, under certain conditions, the difference |r 0 (£ ) - n can be made small 
for £ in fi xed interval. The reflection coeffi cient, K corresponding to q\ tends to zero 
at infi nity. The relationship between m and r\ is studied below. 

The SLR approach should be contrasted with the approach outlined above: Given 
r we approximate it by a smooth magnetization profi le r s . We then construct a smooth 
potential q s with reflection coeffi cient r s . Such a pulse is again implemented by replac- 
ing it with a sum like q t . So there should be no confusion, we denote a sequence of 
boxcar pulses approximating q s hyq sX . Below we show that the difference in the mag- 
netization profile produced by % and mat produced by q sX is controlled by Wqs\-qA\o- 

For our analysis of the hard pulse approximation we consider the one parameter 
family of pulses 

N 

* = J2 ^X[tUA)(/ - J A). € 6 [0, 1], (56) 

with the understanding that q 0 is the sum of 3-pulses given above. It is a simple com- 
putation to show that 



/ 
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is independent of € and 

oo oo 

j ifcCQpA-i y i ?1 (oiv/. 

— oo —oo 

This formula indicates the desirability of using q x to approximate q 0 rattier than a* for 
ane < 1. * 

Let r € denote the reflection coeffi cient defi ned by q. The potential q € is supported 
in the interval (0, NA\. In order to study r € we need to fi nd a formula for _(f ; ATA). 
As is well known, this can be conveniently expressed in terms of a product of 2 x 2- 
matrices. The solution operator at time /<> + A to the ZS-system with potential given 



€ Vx[0, fA )(/-/o) 



is 



( 57 ) 

where 

« = V? 2 + N 2 . (58) 

The solution operator for the ZS-system, at time NA with the potential given by (56) 
is therefore the product of the 2 x 2-matrices: 



It follows easily from (34) that 



^p € (£) -iT*($)1 (59) 



Aft)*-**** 

Xft> • w 

The entries of (57), (59) and (57) are are continuous at € = 0. The functions {A €% B e ) 
are the scattering coeffi ciente for the ZS-system with potential %. If e ^ 0 the asymp- 
toucs established in (65) are applicable. If € =» 0 then A 0 and B 0 can be expressed 
as 

Where P and g are polynomials of degree N - 1. As expected ^ 0 and * 0 are A _1 2jt- 
periodic functions. 

The SLR design process begins with r 0 (£), a rational function of an appropriate 
degree m e l *j, approximating me desired magnetization profi le r (f ) over the interval 
t~ A ~ A l 7t]. The algorithm then determines coeffi cients so that 
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(a) A plot of D\. 




Figure 1 : Plots of the error functions D\ and D 2 for /* e [0, 25]. 

As i with the 1ST approach this problem does not have a unique solution. The different 
solutions are labeled, in essence, by the locations of the zeros of P outside die unit disk 
These correspond to zeros of A 0 in the upper half plane. The questions of principal 
interest are: 

1 . How well does n (f ) approximate r 0 (f ) over the interval [- A"" 1 n 9 A" 1 jr]? 

2. How rapidly does r\ decay outside of [- A" 1 jt, A' 1 jt]? 

3. How are the zeros of A\ in the upper half plane related to those of A 0 7 

The first question admits of a fairly simple analysis. We first consider a sin- 
gle terra in the product expansion (59), exaimning the dependence of the difference 
flC'ji.i.M?) - Oii.o, A<f)!l on (At, A) for £ in the interval [-A-'jt, A~ j tt]. By setting 
£ = Af this reduces to consideration of the differences 



VM 2 + ( 2 



Z*sin|Ml 



(61) 



for f e [-it,n\. As the plots in figures 1 and 2 show, the only way to make both 
differences small is to take ImI small. 

In terms of pulse design, this means that A should be taken to be small which in 
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Figure 2: Plots of the error functions D\ and Z>2 for /z e [0, 0.11. 
turn forces N to be large. A computation with Taylor series shows that 

ih(ti,o = \tiim 2 GQi.t)+( 2 H(t 2 ), (62) 

here F y G t H are entire functions. For a fi xed target magnetization r the f /x* 1 anoearine 
in (59) satisfy ^ * F 6 

|^ / | = 0(A)=a(Ar 1 ). 

For | restricted to si fixed fi nite interval, £ also behaves like OCiV* 1 ). With this restric- 
tion on f both error functions behave like 0(N~ 2 ) as N tends to infi nity. 
We express U\ t & as 

where 

Using this expression and the estimates above one can show that the differences 
Mi(*)-4><t)l. l*itt)-*o«)| 

can be made uniformly 0(N~ l ) for £ belonging to a fixed interval, provided N is 
suffi cientiy large. Note that if the error in the individual terms were of order IT 1 then 
this statement would not be true, as indicated by the well known formula 
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(a) A plot of \rol the 
hard pulse approxima- 



2W" 



(b) A plot of | , an m- 
tcnnediate case. 



(c) A plot of \r\ U the soft 
pulse approximation. 



Figure 3: Plots of the absolute magnetization profi les for various approximations to a 
90°-sinc pulse showing several fundamental periods. 

Thus by taking N large and therefore A small the softened SLR pulse should pro- 
vide a good approximation to the desired magnetization over a fi xed interval where 4) 
does not vanish. For a fixed TV the behavior of ntf) for large |£| is governed by (65). 
Empirically it is well known that r\ is quite small outside of [-A" 1 ^, A" 1 *]. The 
plots in fi gure 3 shows the magnitude of the magnetization profi les for several approxi- 
mations to a 90°-sinc-pulse with a bandwidth of 1000Hz. Note mat r i has considrable 
support outside the fundamental period [-500, 500]. 

The answer to the last question is not immediately clear. According to [10], the 
mmimum energy SLR pulse is attained by placing all the roots of P inside the unit 
circle This accords well with our analysis, as it is zeros of a in the upper half plane, 
which corresponds to the exterior of the unit circle, that increase the energy require- 
ments of an RF-envelope. The asymptotic analysis above shows that zeros of A 0 in the 
upper half plane that are close to the real segment [— A^jt, A" 1 jt] are likely to pro- 
duce zeros of A , in the upper half plane which thereby increase the energy requirement 
of the softened pulse. The behavior of the zeros of A € as * goes from 0 to 1 seems an 
interesting question for further investigation. 

By way of contrast, using the 1ST method one is in direct control of the zeros of 
a m the upper half plane as well as the phase of the magnetization profi le. We close 
this section by stating an estimate which shows that replacing the smooth RF-envelope 
found using the 1ST by a soft pulse approximation q sX does not result in substantial 
errors m the magnetization profi le. 

Theorem 8. Suppose thath\ and h 2 are integrable potentials for the ZS-system, and 
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(63) 
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that ^(f ), i = l,2 solves the ZS-system with potential h it i » 1. 2. Let 
j \h x {s)-h 2 {s)\ds. 6^(0= j \hx(s)\ds % 

-00 -oo 

/oo oo 

1*1 CO - A2(r)k£9. fi + (0 = / \h x (s))ds. 
' 1 

If 

/Ae/i, for f € K f/re differences satisfy the estimates 

(§; 0 - if (£; OO < 2/ o (20/i2BLi)H ± (r)exp(e= fc (O)- (64) 

The proof is in Appendix B. 

In light of these estimates it is clear that by choosing a soft approximation q 8 \ to 
the smooth pulse q s so that Jiq s - q s \ fl^i is small, the reflection coeffi cient r s \ can be 
made uniformly close to r, over any desired interval where a a does not vanish. Once 
again this amounts to using a small value of A in the defi nition of $i . 

9 Appendices 

9.1 Appendix A: Solving the Marchenko equation II 

We now explain the modifi cations to the algorithm presented in section 5 which are 
required if the intermediate potentials {qj) have non-trivial bound states. Suppose that 
<p is a function with two continuous derivatives and -oo < a < fi < oo. In Appendix 
B we establish the following asymptotic formulae and estimates for the scattering data 
(a Ct b c )* for the ZS-system with potential X[a,fi}<P : 

Here T\p,r.A,B] is the piecewise, linear function 

W) = WO *'-^ 0 '- 0 . (66) 

These formulae hold uniformly in the upper half plane with the implied constants in the 
0-terms bounded in terms of the %\[a t £])«norm of q>. The formula for a e shows that 
the zeros of this function, in the upper half plane are constrained to lie in a fixed ball 



(65) 
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of radius Ro. Dividing these formulas gives an asymptotic formula and estimate for the 
reflection coeffl cient, r c : 

e*Wr c {$) = ^f M *to*wm) + O (^^2) - (67) 

The error term in again uniformly estimated in the upper half plane outside the ball of 
radius R 0 + 1 (say). Observe that both r c and d^r c belong to L 2 (JK). Thus the results 
of Beals and Coifinan show that the inverse scattering transform can be applied to this 
data. 

We can now complete the proof of Theorem (7) 

Proof of Theorem (7). if q has support in (-oo, /+] then the scattering coeffi cient b 
has an analytic extension to the upper half plane. This is because = [e~'*' t 0]t 
for t > /+ and therefore equation (22) shows that + 

^) = ^2i-tt;^)e- |?r+ . 

If r has a meromorphic extension to the upper half plane then the integrand in the 
Marchenko equation, defi ned in (36), has a different representation. If the zeros of a in 
the upper half plane have imaginary parts less than no then 

—00 

see [1]. Provided that the potential has support in a positive half line, there is an anal- 
ogous formula expressing g. the kernel function in the left Marchenko, as an integral 
overthe line {f - /n 0 : 170 e K} : 

00 

2* J a(?-/i, 0 ) *" (69) 

— OO 

Using the asymptotic formula given in (65) and the asymptotics for # 21 _ that fol- 
low from (16) and (18), we conclude that e*S'+r satisfies the hypotheses ofthe Paley- 
Wiener Theorem II in the half plane Im£ > » 0 (see Appendix B). This function is 
therefore the Fourier transform of a function with support in (-oo, 0]. Hence / is sup- 
ported in (-00, 2£j. This argument applies mutatis mutandis to the kernel ofthe left 



Remark 9. If a has only simple zeros then (68) implies that the norming constants are 
given by 

There is a similar formula for the [C, J. Given that the potential has support in an appro- 
priate half line, the formula* (68) and (69) defi ne the kernel functions for the Marchenko 
equation whether or not the zeros of a are simple. 
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Since the potentials [qj] have bounded support and q is assumed to have at least 
two derivatives the estimates in (65) apply for each j > 1 . The scattering coeffi cients 
{(fly. bj)) are entire functions. The asymtptotic formula for aj implies that is has at 
most finitely zeros in the upper half plane. Denote these zeros by :i = \,... 9 Nj) 
At step j we have determined 

*7 = X<-^,,<T">r* + X [t O-i> %ooy q. 
Equation (39) implies that 

OO OO 

j | 9 (/)| 2 < // = i f log(l+|r(£)| 2 )</|. 

— OO — oo 

It is not diffi cult to show that (39) can be extended to potentials of the form 
with q twice differentiable, so that 

oo oo tfj 

j \qj(t)\*dt = ± j log(l + |ry(f)| 2 )^ + 4X;imf / (/) . 

-oo -oo '= 1 

Using these formula? gives 

This relation gives an upper bound on the imaginary part of any root of ay : 



-00 \ J / -oo 

oo 



(71) 



We fi x tft to be a non-negative number, larger than max, j {Im f/^J or zero if this max- 
imum is zero. 

In light of the upper bound (71) on the imaginary parts of the zeros of the functions 
{fly}, we can use the estimates in (65) to further conclude that the integrands {/y : j > 
1} are uniformly bounded in t2fff~ ! \ ^"^l. Indeed each is the sum of a piecewise 
linear function of the type defi ned in (66) and a continuous function with support in 
(-co, 2t^"^] which vanishes at -co. In light of Theorem 5 this establishes the claim, 
made in section 5, that mere is a A > 0 so that 
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These observations show that the algorithm in section 5 applies to the case where 
the intermediate potentials have bound states so long as we use (68) ((69)) to compute 
*e {/)} ({g>)). This formula requires the analytic continuation of the scattering data 
into the upper half plane. Since we know, a priori that a, and bj have such extensions, 
the estimates in (65) apply to show mat the Cauchy integral formula gives: 



50 (M?> -D<*? 



— OO 

(72) 

2*' J f-(f + *jj) 
— oo 

To obtain an algorithm that succeeds in solving the Marchenko equations in all 
cases we need only modify Step 1 in the earlier algorithm: 

Step 1 • : Use (28) to determine if the reflection coeffi cient for q, has poles in the up- 
. per half plane. If so extend (a,, bj) to the line Imf = 170 and use (68) ((69)) 
to determine the kernel runction(s) fj (gj) supported in the ray (-oo, 2f^~ 1> ] 

o^sr^oo)). 

We close this appendix with some remarks on the implementation of this algorithm. 

1 . It is well known mat it is not necessary to solve the 2 x 2-matrix system (37) 
to fi nd 3$2/ • By substitution this can be reduced to solving a scalar integral equa- 
tion, see [1]. 

2. Given that the initial data r (or m°°) is assumed to be smooth, with bounded 
support there is no problem, in principle, in determining q^. The scattering data 
for these potentials can be determined by numerically solving the ZS-system 
with the given potentials. The method used in section 8 gives alternate algebraic 
approach to the approximate determination of the scattering coeffi cients. 

3. In the subsequent steps of the algorithm we are working with potentials with 
bounded support but jump discontinuities. The difficulties that this might en- 
gender, like the Gibbs phenomenon, can be avoided by using the asymptotic 
expansion (67) in (68). That is, we set 

The contribution of ^-"^-"K^^O-o^^^O-Ojj to /} can be exactly deter- 
mined. It is a function with support in [2^"°, 2/£ _I) ]. The function r, is 
0(11 1~ 2 ) as HI goes to infinity. This estimate holds uniformly in strips with 
bounded imaginary part in the upper half space. Thus in the computation of the 
approximate Fourier inverses 



K 



231 

-R 
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we can avoid bothersome Gibbs-type artifacts. If q is smooth then there is a com- 
plete asymptotic expansion for each r, which only depends on the (in principle 
known) values of the derivatives [q W Using further terms in this expan- 

sion can produce even more rapidly convergent correction terms. This should be 
done with care as it requires the computation of derivatives of q which could in 
turn lead to numerical instability. 

4. After Step 0 we are looking for spatially limited potentials. The input data for 
Step 1 is essentially in the frequency domain. Since we know t^] we can 
choose the sample spacing in the frequency domain so as to avoid aliasing errors 
in the construction of these potentials. 

5. The Cauchy integral in (72) is a shift invariant filter with transfer function 
2^^/x<-oo t 0](0e /,, . As such it is easily implemented in practice. 

6. After a potential q with reflection coeffi cient r and no bound states is found, 
it is a relatively simple matter to add in bound states. An algorithm to do this 
is described in section 6 of [3]. One possible usage for bound states is to re- 
duce the effective support of the potential. To that end one might chose points 
{£i • . . ■ . %n) in the upper half plane and norming constants {Cj } so that sum 

7=1 

approximates 

oo 

• m = h J r ^' d ^ 

—00 

. over some ray [7\ oo). With these choices the effective support of / defined 
in (36) can be made smaller than that of r and thereby a potential with smaller 
(effective) support,' producing the same magnetization profi le can be found. This 
is of course at the expense of 

N 

4£lm* y 

-M 

in additional energy. In section 6 we saw that the effective support of q de- 
termines the rephasing which is needed to obtain the specifi ed magentization 
profile. 

7. While the bound states do not affect the observable magnetization they do' 
affect the spins themselves. Bound states might prove useful in designing pulse 
sequences intended to manipulate higher order quantum coherences. 



8. The computation in (28) has an integer valued result Using the asymptotic 
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formula (65) it is not diffi cult to give an estimate for it so that the fi nite integral 

R 



«>,(£)• (2*0 



/ 



suffi ces to determine N, . 



9.2 Appendix B: Mathematical proofs 

The Paley-Wiener Theorems: The "Paley-Wiener" Theorem is used several times 
in this paper. This is really a collection of results relating the support properties of a 
function defi ned on R to the analytic continuation of its Fourier transform. A complete 
discussion can be found in [17]. Two results adequate for our applications are: 

Theorem 9 (Paley-Wiener Theorem I). A Junction f in L 2 (M) has bounded support 

in the interval [— L , L] if and only if its Fourier transform has an analytic extension to h 

the entire complex plane satisfying: 

l/(!+in)|<Ce^, (73) 

for a constant C. 
and 

Theorem 10 (Paley-Wiener Theorem II). A Junction f in £ 2 (R) has support in the 
ray (—oo, 0] if and only if its Fourier transform has an analytic extension to the upper 
half plane satisfying 



oo 

j \?(M+m 2 <Cand 

— oo 
oo 

Urn j |/«+/i?)-/tt)I 2 = 0. 



(74) 



The proof of Theorem 5: Assuming that //(2t) < 1 We need to show that the 
Marchenko equation, 

00 

W f (5) - fw t ts + x)&(t + x)dx = -®(s + /), (75) 
i 

can be solved by the iteration defined in (44) and the solution satisfies the estimate 
in (45). The proof is a simple induction argument Clearly ||3£?(y)| < M/(t 4- s). 
Assume that M 

fr)|<M/(» + s)[l + //(20 + .../^(2f)]. (76) 



-~ £—1. o i — -r v « o .rti — n TT/ 5 -— i 
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It follows from this assumption that 

oo 

|3t/ +l < M f (f + 5) + [l + //(2f) '+ . . . /J(2<)] x J M/(t +x)W(t + *>!!<**■ 

i 

The induction hypothesis (76), with y replaced by J + 1 follows easily from this esti- 
mate. Using a similar argument we can show that, for j > 1 

-3tf-*(s)| <MfV+s)If{20J (77) 

If J/(2r) < 1 then this estimate shows that {9f/ ) is a uniformly convergent sequence 
for s € (/, oo). The limit, 3tJ° is evidently a solution to (75). We can also pass to the 
limit in (76) to conclude that 3C f° satisfi es (45). 

The proof of (65): We use the integral equations (16) defl ning the components of 
*,_.Let 

/iG;0 = *u-(*;Oe' f '. 

It follows from (15), (16) and (19) that 

} 

a e (g ) = 1 + / (5; oo. s)/i (5; *)<&. 

(78) 

a 

In light of the simple form assumed by the potential we see, integrating by parts, that 

= y«xwi« I 2J1 / — 271 — dx J 

(79) 

Since (x — s) > 0, in the integral, this expression remains uniformly bounded for £ in 
the upper half plane and implies that 

M«; oo.,) = -V* W*> [ ^g>-^> + o (-J-p)] . (80) 

This formula holds for £ with non-negative real part with the implied constant in the 
0-term only depending on the <@ 2 [cx, 0]-norm of y>. 
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Using this formula for M\ (£; co, s) in (78) gives 




(81) 



This completes the analysis of a c \ for b c we integrate by parts again to get 




+ °Griip)- 



(82) 



The proof is completed by observing that 

~ Be~ l * - Ar'#* B-A e - ,w * - 




If y has k + 1 derivatives then we can integrate by parts in (79) A>times to get an 
asymptotic expansion for M\ with an error term of order 0(|f |~ ( * +I) ). This can in 
turn be used in (16) and (78) to obtain higher order asymptotic expansions for a c and 
b c . For smooth potentials, vanishing at infi nity, expansions of this type are well known 
in the inverse scattering literature. 

The proof of Theorem 8: The proof of mis result is a small modifi cation of the proof 
of the estimate (18) given in [1]. Let a denote the 2 x 2-matrix 



We now restrict our attention to the — case, the + case is handled by an essentially 
identical argument 




and let 



For real f and / the matrices e** ttia are unitary and therefore 
The vector valued functions satisfy the equations 
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The differential equation in (83) is equivalent to the integral equation 

-00 

/ 

J e L-<AJ-AS><*> 0 \ F * M*' 

— OO 

This is a Volterra equation which can be solved by the following iteration: 
-oo 

I 



for y > 1. 

First observe that the operator norm defi ned by the Euclidean norm satisfi 

0 al 
-or* 0 



= M 



and therefore (1 8) implies that 

Wo Oil < 2/ 0 (OA 2 || L i)/r-(0; 

we make the inductive assumption that 

il^JG; Ofl < 2/o(p 2 IUO"-(0 (g ^ 0) * j • 
From these assumptions and(85) it follows that 

Bi>7 +1 «;OD<2/o(DA2llxO + | |Aifr)|/r(5) Pfc (g ~f )} * 

\ -oo, [k=0 

The defi nition of Q~ implies mat 

|Ai(s)| = a 4 fi-(s) 

and therefore (87) implies that 

IIDJ+.tf ; Ofl < 2/ 0 (||A 2 ii / ,,)h-(o |^i + j a s ^C ^^y ^ <fcj 



•q» iU- o »-?. t t 1 : *• * . er, „a Br. 
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thus completing the proof of the induction step. As 

lim D~ = £>", 

letting 7 tend to infi nity completes the proof of the theorem. 
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I CLAIM: 

1 . A method of synthesizing minimum energy pulses for a given magnetization 
profile in a magnetic resonance imaging system, comprising the steps of: 

determining minimum energy for an RF envelope producing the given 
magnetization profile; 

using an inverse scattering transform to determine a minimum energy smooth 
pulse RF envelope that approximates the energy of said minimum energy RF 
envelope; and 

generating pulses with energy that approximates said minimum energy smooth 
pulse RF envelope. 

2. A method as in claim 1, comprising the step of determining auxiliary 
parameters including bound states and norming constants for reduced scattering data 
of a potential for the minimum energy used in said inverse scattering transform. 

3 . A method as in claim 1 , wherein the energy required for the minimum energy 
smooth pulse RF envelope is independent of the phase of transverse magnetization. 

4. A method as in claim 1, comprising the additional step of increasing the 
energy of the smooth pulse RF envelope in order to reduce rephrasing time. 

5. A method as in claim 2, wherein the step of determining the minimum energy 
smooth pulse RF envelope comprises the step of solving left and right Marchenko 
equations by iteration on maximal intervals to determine the reduced scattering data 
and iteratively removing parts of the reduced scattering data produced by already 
known parts of the potential. 

6. A method as in claim 1, comprising the step of replacing the smooth pulse RF 
envelope with a soft pulse approximation to said smooth pulse RF envelope. 

7. A method of generating a desired frequency dependent excitation in a 
magnetic resonance imaging system using minimum energy selective RF pulses for a 
given magnetization profile corresponding to said desired frequency dependent 



excitation, comprising the steps of: 

determining minimum energy for an RF envelope producing the given 
magnetization profile; 

using an inverse scattering transform to determine a minimum energy smooth 
pulse RF envelope that approximates the energy of said minimum energy RF 
envelope; 

generating pulses with energy that approximates said minimum energy smooth 
pulse RF envelope; and. 

applying the pulses to said magnetic resonance imaging system to obtain a 
desired selective excitation. 
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